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Observation of (3686) — e*e^x,; and Xes > ete J/w 
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Using 4.479 x 10° 4(3686) events collected with the BESIII detector, we search for the decays 
(3686) > e*e xc; and xe; — ete J/p, where J = 0,1,2. The decays (3686) > e*e xey and 
XcJ — €*e- J/w are observed for the first time. The measured branching fractions are B(w(3686) > 
ete xej) = (11.7 2.53: 1.0) x 1074, (8.62:0.32: 0.6) x 1074, (6.93: 0.5 0.6) x 1074 for J = 0,1,2, 
and B(xc; + ete J/w) = (1.51 + 0.30 £0.13) x 1074, (3.73 + 0.09 + 0.25) x 107%, (2.48 + 0.08 + 
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B(V(3686)—e* e^ xes) 


0.16) x 107? for J — 0,1,2, respectively. The ratios of the branching fractions a ICEL) = 3x.) 


B(xoy et e J/ip) 
B(xej—7J/w) 
e pair are determined for 7(3686) > e 


and 


+ E 


€ 


PACS numbers: 13.20.Gd, 13.40.Hq, 14.40.Pq 


Study of electromagnetic (EM) Dalitz decays i, in 
which a virtual photon is internally converted into an 
e*e- pair, plays an important role in revealing the struc- 
ture of hadrons and the interactions between photons 
and hadrons [3l Such decays are widely observed in the 
light-quark meson sector, for example, n! — ye*e- ,9' > 
we*e-, and ó — nete~ However, the analogous 
transitions in charmonium decays have not yet been stud- 
ied. Although the potential quark model has successful- 
ly described the low-lying charmonium states with high 
precisions, there are still puzzling discrepancies in the 
decay branching fractions B(~(3686) — YXeJ) between 
the experimental results |3] where the higher-order mul- 
tipole amplitudes are ignored and the various theoretical 
predictions . Throughout this Letter, y. 7 refers to 
Xc0,1,2. While recently the BESIII experiment confirms 
that the contributions from the higher-order multipole 
amplitudes in (3686) — yxez are small ER the E1 con- 
tribution is dominant. Therefore, it is of great interest 
to measure the EM transition 7(3686) — ete~x-7 and 
Xer > ete” J/y. 

The EM Dalitz decays in charmonium transitions, such 
as (3686) — e*e^x«j or Xez > e*e J/, have ac- 
cess to the EM transition form factors (TFFs) of these 
charmonium states. The q?-dependence of charmonium 
TFFs can provide additional information on the inter- 
actions between the charmonium states and the elec- 
tromagnetic field, where q? is the square of the invari- 
ant mass of the e*e^ pair, and serve as a sensitive 
probe to their internal structures. Furthermore, the q?- 
dependent TFF can possibly distinguish the transition 
mechanisms based on the cc scenario and other solutions 
which alter the simple quark model picture. We em- 
phasize that the q?-dependent TFF can also serve as an 
useful probe for exotic hadron structures based on differ- 
ent models. One example is that with the precise mea- 
surement of the radiative decay of X (3872) > e*e- J/w 
and X(3872) — e*e w(3686) in the future, we can 
pin down the intrinsic structure of X (3872) by compar- 
ing the experimental measurement of the q?-dependence 
of TFF with different model calculations. The nature 
of X(3872), namely whether it is a compact charmoni- 
um, multiquark state with quark clustering, or hadronic 
molecule ire can possibly be disentangled by the q?- 
dependence of its TFF. 

In this Letter, we report the observation of the 
EM Dalitz decays (3686) — e*e^x,; and Xe — 
ete J/wW by analysing the cascade decays (3686) — 
ete Xej,XcJ — ^*J[v and v(3686) > YxXeI,Xer > 


are also reported. Also, the o values of helicity angular distributions of the 


€ Xe1,2 and xa; + ete J/v. 


e*e- J/ip, respectively. Here, the J/w is reconstructed 
in its decay to an ete~ or ppt pair. The two cascade 
decays studied have the same final state: four leptons 
and a single photon. The analysis uses a data sample of 
4.479 x 108 (3686) events (14, [15] taken at a center-of- 
mass energy y/s = 3.686 GeV collected with the BESIII 
detector lid operating at the BEPCII storage ring in 
2009 and 2012. In addition, a data sample corresponding 
to an integrated luminosity of 44 pb~!, taken at a center- 
of-mass energy ys = 3.65 GeV , is used to estimate 
the background from continuum processes. 

The BESIII detector has a geometrical accep- 
tance of 93% of the total 4r solid angle. A small-cell 
helium-based main drift chamber (MDC) provides mo- 
mentum measurements of charged particles with resolu- 
tion of 0.5% at 1 GeV/c. The MDC also supplies an 
energy loss (dE/dx) measurement with a resolution bet- 
ter than 696 for electrons from Bhabha scattering. The 
time-of-flight system (TOF) is composed of plastic scin- 
tillators with a time resolution of 80 (110) ps in the bar- 
rel (endcaps) and is used for charged particle identifi- 
cation. The CsI(TI) electromagnetic calorimeter (EMC) 
measures 1 GeV energy photons with a resolution of 2.596 
(596) in the barrel (endcaps) region. 

Monte Carlo (MC) simulations are used to estimate the 
reconstruction efficiencies and study the backgrounds. 
The signal MC samples are generated using EVTGEN 
using a q?-dependent decay amplitude based on the as- 
sumption of a point-like meson, as described in Ref. Rd, 
and an angular distribution based on that observed in 
data. An MC sample of generic 7)(3686) decays, the so 
called “inclusive MC sample”, is used for the background 
studies. The production of the 7(3686) state is simulat- 
ed by the KKMC generator. The known decay modes 
of the 7)(3686) are simulated by EVTGEN according 
to the branching fractions reported in PDG [3], while the 
unknown modes are simulated using the LUNDCHARM 
model. 

Each charged track is required to have a point of clos- 
est approach to the interaction point (IP) that is less 
than 1 cm in the radial direction and less than 10 cm 
along the beam direction. The polar angle 0 of the 
tracks must be within the fiducial volume of the MDC 
(|cos0| < 0.93). Photons are reconstructed from isolat- 
ed showers in the EMC which are at least 20? away from 
the nearest charged track. The photon energy is required 
to be at least 25 MeV in the barrel region (| cos 6| < 0.8) 
or 50 MeV in the endcap region (0.86 < | cos0| < 0.92). 
In order to suppress electronic noise and energy deposi- 
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tions unrelated to the event, the time after the collision 
at which the photon is recorded in the EMC must be less 
than 700 ns. 

Candidate events are required to have four charged 
tracks, with a sum of charges equal to zero, and at least 
one photon. The tracks with momentum larger than 
1 GeV/c are assumed to be leptons from J/w decay. 
Otherwise they are considered as electrons from the v 
or X-y decay. Leptons from the J/wv decay with EMC 
energy larger than 0.8 GeV are identified as electrons, 
otherwise as muons. The J/* signal is identified by re- 
quiring the invariant mass of the lepton pair to be in the 
interval [3.08, 3.12] GeV/c?. A vertex fit is performed 
on the four charged tracks to ensure the tracks origi- 
nated from the IP. In order to reduce the background 
and improve the mass resolution, a four-constraint (4C) 
kinematic fit is performed by constraining the total four 
momentum to that of the initial beams. If there is more 
than one photon candidate in an event, all the photons 
are individually fit with the four leptons in the kinematic 
fit and only those with a fit x? < 40 are retained. If two 
or more photons pass this criterion, only the one with 
the least x? is retained for further analysis. 

A study of the (3686) inclusive MC sample shows 
that, after applying the above selection criteria, the main 
background comes from (3686) > YxXer,Xer > vyJ/v 
decays, where one photon converts into an ete~ pair in 
the detector material. To suppress this background, a 
photon-conversion finder is applied to reconstruct 
the photon-conversion vertex. The distance from the 
point of the reconstructed conversion vertex to the z 
axis, R,y, is used to distinguish the photon conversion 
background from signal. By studying the MC sam- 
ples (3686) > «xc, xcu > YJ/v, the peaks around 
Rey = 3 cm and Rzy = 6 cm match the positions of the 
beam pipe and the inner wall of the MDC , respec- 
tively. We remove the events in 1.5 cm< Ryy < 7.5 cm 
to suppress the y conversion background. With this re- 
quirement, the y conversion background is negligible for 
the decays (3686) + e*e^ xez and is at the few percent 
level for the decays yey > e*e- J/v. 

To remove the backgrounds from decays (3686) — 
1/19 Jb, n/n? — yete, which have the same final state 
as signal events, a requirement 0.16 < M(yete—) < 0.50 
GeV/c? is applied. By studying the data collected at 
Vs = 3.65 GeV, the contribution from the continuum 
process is found to be negligible. 

Figure [1] shows the scatter plot of M(»yJ/wv) versus 
M(e*e-J/v) for the selected events from data; the 
corresponding one-dimensional projections are shown in 
Fig. Clear x.y signals are observed in the M(yJ/w) 
and M(e*e- J/v») distributions, corresponding to the de- 
cays 1/(3686) > e*e^ xez and Xez — e*e- J/y), respec- 
tively. The study of (3686) inclusive MC samples in- 
dicates that the dominant background is from the decay 
(3686) > vm JI, JI > (ygsn)U L^, where yfsr is 


a photon due to final-state radiation; these events accu- 
mulate at M(ete~ J/w) ~ 3.6 GeV/c?. 


3.67 


M(y J/v) (GeV/c?) 


3.2 34 36 
M(e'eJ/y) (GeV/c?) 


FIG. 1. (color online) Scatter plot of M(yJ/w) versus 
M(e*e- J/w) for data. The horizontal red dashed lines and 
vertical blue dashed lines indicate the positions of the Xes 
masses in the M(yJ/w) and M(e*e- J/w) distributions, re- 
spectively. 
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FIG. 2. (color online) Data (points with error bars) distri- 


butions of (left) M(yJ/w) and (right) M(e*e- J/w). The red 
solid curve is the overall fit result, the green long-dashed curve 
is for the background (left) 4(3686) — YXco, Xeo > ete J/w 
and (right) (3686) — ete xco,xco — yJ/w, the blue 
dashed curve is for QED background, and the pink dashed- 
dotted curve in right plot is for the backgrounds from w(3686) 
decays. 


Separate unbinned maximum likelihood fits are per- 
formed on the M(yJ/) and M(e*e- J/1) distributions 
to extract the signal yields. We use the signal MC- 
determined shape, convoluted with a common Gaussian 
function, to describe the shapes of y-y signals. The 
Gaussian function parametrizes any resolution difference 
between the data and MC simulation and its parameters 
are determined from the fit. 

Two background components are considered in the 
fit to the M(yJ/p) distribution. The first background 
is from the decay (3686) > YXco, Xeo — ete J/y, 
which corresponds to the peak at the lower edge of the 
M(yJ/w) region; it is described by a MC-determined 
shape with a fixed number of events based on the branch- 
ing fraction obtained in this analysis. The second one is 
related to QED background (ete~ — £*£ ,£ = e,u,T) 
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and is described by a first-order polynomial function in 
the fit. 

In the fit to the M(ete~ J/w) distribution, three back- 
ground components are considered. The first two are 
from the decay (3686) — e* e^ xco, Xco — *J/A, which 
corresponds to the enhancement at the lower edge of 
the M(ete~ J/w) fit interval, and QED processes; the 
way these components are dealt with in this fit is anal- 
ogous to the way they are dealt with in the M(4J/w) 
fit. The third background component is from inclu- 
sive 1/(3686) decay, which includes the dominant one of 
(3686) — ntr J/v,Jlb — (versn) l7. decays and 
a small fraction from (3686) > oixej,Xej > voJ/ v, 
where y> converts into an e*e^ pair. In the fit, the 
shape of the third background component is assumed to 
be that reconstructed in the inclusive MC sample with 
the normalization determined from data. The fit results 
are shown in Fig. 2] and the corresponding signal yields 
are summarized in Table [I| For the six observed decay 
modes, the statistical significance of the yields are all 
larger than five standard deviations. 

The branching fractions B((3686) — e*e^x.7) and 
B(xc; > ete” J/)) are calculated according to 


Nsig 
Na(a686) * € * Bradiative * B(J/ > L07)' 


where Nsig is the corresponding number of signal events 
extracted from the fit, Ny36g6) is the total number of 
(3686) events, e is the selection efficiency determined 
from the signal MC samples, Byadiative is the branching 
fraction of the radiative transitions (3686) > YXeJ or 
Xer > *J[b, and B(J/ — IIT) is the decay branching 
fraction of J/w — l*1^. All the branching fractions used 
are taken from Ref. |3]. The resultant branching fractions 
of (3686) — e*e^x,; and Xez — ete J/w are listed 
in Table [I] 

Figure [3| shows comparisons of the q distributions 
in data and MC simulation for the decays (3686) > 
ete Xe1,2 and xa, > ete J/w, where the xc and xc» 
signals are extracted requiring a mass within [3.49,3.53] 
and [3.54,3.58] GeV/c”, respectively; with these criteria 
the backgrounds are expected to be less than 2%. The 
data are in reasonable agreement with the MC simulation 
generated using the model described in Ref. Ro. 

The systematic uncertainties for the branching frac- 
tion measurement arise from the following sources: track 
reconstruction, photon detection, kinematic fitting, J/1 
mass criteria, M (ye* e^) requirement, y conversion veto- 
ing, fit procedure, angular distributions, the total num- 
ber of (3686) events and the branching fractions of the 
cascade decays. All uncertainties are discussed in detail 
below. 

'The difference in the tracking efficiency between data 
and the MC simulation, for each charged track, is esti- 
mated to be 1.096 ld, which results in a 4.0% system- 
atic uncertainty for all modes. The uncertainty on the 
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FIG. 3. Data to MC simulation comparisons of q distribu- 
tion for the decays (a) (3686) — e*e xa, (b) w(3686) > 
e*e xe, (c) Xa > ete J/y and (d) xe > ete J/w. The 
points with error bars are data and the red histograms are for 
the signal MC simulation. 


photon-detection efficiency is derived from a control sam- 
ple of J/ — p?r? decays and is 1.0% per photon [25]. 

In the 4C kinematic fit, the helix parameters of charged 
tracks are corrected to reduce the discrepancy between 
data and the MC simulation as described in Ref. [26]. 
The correction factors are obtained by studying a control 
sample of (3686) > ntr J/w, J/p — Il decays. To 
determine the systematic uncertainty from this source, 
we determine the efficiencies from the MC samples with- 
out the helix correction; the resulting differences with 
respect to the nominal values are taken as the systemat- 
ic uncertainties. 

The uncertainty associated with the J/wW mass require- 
ment is 1.096, which is determined by studying a control 
sample of 1»(3686) > nJ/V,r — yy (where one y under- 
goes conversion to an ete~ pair) or n + yete~ decays. 
The systematic uncertainty related to the M(ye*e7 ) in- 
terval used is studied by varying the edges of the interval 
by +5 MeV/c?. The largest difference with the nominal 
value is taken as the systematic uncertainty from this 
source. 


To study the systematic uncertainty related to the y 
conversion background veto, we compare the efficiencies 
of y conversion veto between data and the MC simu- 
lation in control samples of (3686) — YxXe1,2; Xe1,2 — 
ete J/w decays. The efficiency of the y conversion veto 
is the ratio of the signal yields determined by fitting the 
M (e*e-) distribution with and without the y conversion 
veto applied. A relative difference between data and sim- 
ulation of 1.496 is found and assigned as the systematic 
uncertainty. 

'The sources of uncertainty in the fit procedure include 
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TABLE I. Signal yields, detection efficiencies, the branching fractions and the ratios of the branching fractions. Here the first 


uncertainty is statistical and the second systematic. 


B(v(3686)3e" € xey) B(xcjy—e* e J/v) 


IU RuaRansste-.e .S — Rt. .sta— F/4hN  — 
Mode Yields  Efficiency(96) Branching fraction CIEE ESSE B TAM 
(3086) > ee xao 48 10 6.06 üL7TEX25110)x10-7 (9.44 1.9 -E 0.6) x 107 = 


(3686) — e*e^x« 873 + 30 5.61 (8.6 x: 0.3 
v(3686) + ete-xe3 227416 3.19 (6.9 3: 0.5 
Xoo > ete- Jf 56411 6.95 (1.51 + 0.30 


Xa  ete- Jf 
Xc2 > ete Jib 


1969 + 46 10.35 (3.73 + 0.09 
1354 + 39 11.23 (2.48 + 0.08 


the fit range and the signal and background parametriza- 
tion. The uncertainty related with the fit range is ob- 
tained by varying the limits of the fit range by +5 
MeV /c?. The largest difference in the signal yields with 
respect to the nominal values is taken as the systematic 
uncertainty. In the nominal fit, the signal shapes are de- 
scribed with the signal MC simulated shapes convoluted 
with a Gaussian function. An alternative fit is performed 
by fixing the signal shapes to those of MC simulation. 
The resultant change in the signal yields is taken as the 
systematic uncertainty. The uncertainty associated with 
the background shape is estimated by an alternative fit 
replacing the first order polynomial function with a sec- 
ond order polynomial function for the background shape, 
the resultant change in the signal yields is taken as the 
systematic uncertainty. 


The distribution of ete~ pair's helicity angle in its 
mother rest frame e+e- may affect the detector efficien- 
cy, where e+e- is the polar angle of e*e^ pair in the 
colliding beams rest frame with the z axis pointing in 
the positron beam direction. The efficiency corrected 
cos e+e- distributions are shown in Fig. [4] for the de- 
cays (3686) > e*e^ xc» and Xe1,2 + e*e-J/wv; each 
distribution is fit with a 1 + a cos? 0,.,- function. The 
resultant o values are —0.6+0.2, —0.9+0.3, 0.00.2 and 
0.5 + 0.2 for the decays (3686) + e*e^ xa, v(3686) > 
ete xe, Xa > ete J/y and Xe2 > e*e- J/yY, respec- 
tively. The measured a central values are incorporated 
in the nominal MC simulations. To take into account 
any effect on the detection efficiencies due to an incor- 
rect simulation of the cos @,+,- distribution, alternative 
MC samples are generated with a varied by +1 stan- 
dard deviation and the efficiencies are determined. The 
differences with the nominal efficiencies are taken as the 
systematic uncertainties from this source. In the decays 
(3686) 4 e* e^ x«o and xeo — ete J/w, the cos Oe+¢- 
distribution is not extracted directly from the data due to 
the limited statistics. The theoretical expectations for a 
are 1 and 0 for (3686) — et e~ xXco and Xeo > ete J/wv, 
respectively, which are used to generate the nominal MC 
simulation. The systematic uncertainty is estimated us- 
ing the difference in efficiency when alternative MC sam- 
ples with a = 0 for 4/(3686) > e*e^ xco and a = 1 for 
Xco > ete J/w are used. 

The total number of (3686) events is measured to 
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FIG. 4. Distributions of efficiency corrected cos0,...— for the 
decays (a) (3686) + e* e^ xa, (b) v(3686) > e*e^ xez, (c) 
Xa ete Jh and (d) xca > ete Jfu. The red line is the 
fit to 1 + a cos? 0,.,.. 


within 0.7% by using the inclusive hadronic events [14 


[15]. The uncertainties of the branching fractions in the 


cascade decays are taken from Ref. [3. 


The effect of other potential systematic uncertainty 
sources are considered, such as uncertainties on the gen- 
erated q distributions, the trigger efficiency, and the sim- 
ulation of the event time, but are all found to be negligi- 
ble. Table[IT] summarizes all individual systematic uncer- 
tainties, and the overall uncertainties are the quadrature 
sums of the individual ones, assuming they are indepen- 
dent. 


In summary, using a data sample of 4.479 x 105 (3686) 
events collected with the BESIII detector operating at 
the BEPCII collider, the decays (3686) > e*e xe 
and x«j — ete J/w are observed for the first time, and 
the corresponding branching fractions are measured and 
the values are given in Table [I] The ratios of branching 


B(1»(3686)—e* e^ xes) B(xer vet e7 Jfi) 
BGS) 7x) and =a a are also 
obtained by incorporating the BESIII results of the prod- 


uct of branching fractions B(v(3686) > yxc;): B(xez > 
yJ/w) in Ref. [8], as listed in Table [I| The common sys- 
tematic uncertainties related to efficiency and branch- 
ing fractions cancel in the calculation. The measured q? 


fractions 


201711.02286v1 


chinaXiv 


TABLE II. Summary of systematic uncertainties (in 96). 
(3686) 2 e*e xej  Xegz > ete J/v 


Tracking 4.0 4.0 4.0 4.0 4.0 4.0 
Photon 1.0 1.0 1.0 1.0 1.0 1.0 
Kinematic fit 1.6 1.4 1.4 1.8 2.2 2.4 
J/w mass window 1.0 1.0 1.0 1.0 1.0 1.0 
M(yete-) 2.7 1.2 1.0 0.7 2.2 0.4 

^y conversion vetoing 1.4 1.4 1.4 1.4 1.4 1.4 
Fit Range 2.2 0.2 0.3 4.7 0.1 0.2 
Signal shape 0.4 0.1 0.1 2.2 0.2 0.5 


Background shape 2.2 0.2 0.3 0.1 0.1 0.2 
Angular distribution 3.9 2.1 3.3 3.6 1.6 1.0 
Number of (3686) 0.7 0.7 0.7 0.7 0.7 0.7 
Branching fractions 4.8 3.6 5.5 2.8 3.3 3.5 

sum 8.9 6.5 8.1 8.5 6.6 6.3 


distributions are consistent with those of the signal MC 
simulation based on the assumption of a point-like meson 
2d. This first observation of the q?-dependent charmo- 
nium EM Dalitz transitions can help understand the dis- 
crepancy between the experimental measurements [3] and 
the theoretical predictions Pin of the (3686) + yxe7 
branching fractions. The experimental methods applied 
here for the first study of charmonium Dalitz decays are 
likely to be of use for similar studies of the X (3872). It 
is hoped that this experimental work will spur new theo- 
retical development on use of charmonium Dalitz decays 
to address questions such as the nature of exotic charmo- 
nium. 
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